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Ilea-.    Low-li       nperat    ai       profile   were  con - 
itructe   for   ',  cas   of  free  convection  using  t'ae  ;:ata 
for  O'Nei   ,     aska,  during  July  a:   Augi  it  L956.   Basec 
pon  tl   ■  •>.  '  ision  for  the  uormalized  logarithmic  wind  shear 
first   ■    t<   bj      on  and  later  refine   j  Panofsky,  a 
theoretica   formula  for  tl     aeon  profile  number  as  a  func- 
tion of  the  lichardson  it :r  ''a.;  aerive*  ,  and  /alnes  of 

the  Deacon  profile  number  n/ere  computed.   One  of  the  para- 
meters Altering  into  til   theoretical  formula  is  the  ratio 
of  the  e   .  i  Lff  livities  for  heat  and  momentum.   This  para- 
meter was,  i  ■  t  rn,  computed  from  Priestley'   expression 
for  the   Lm  nsionless  neat  flux  for  free-convective  cases. 
In  using  o  lerved   '  t<   ata  from  the  mean  profile  in  or  er 
to  verif   tl   t'n  oretical  computations  of  _  '  ,  some  marked 
discrepancie;  ccc  rre<  a  ove  the  100  cm  level.   These  were 
due  to  Inconsistent  ;ind  speed  readings,  and  It  ,;as  necessary 
to  employ  control  ;:ata  based  on  neutral  profiles  to  correct 
the   in   pe<   .   When  this  was  done,  the  theoretical  and 
o  servi   >eacon  prpfile  numbers  were  Ln  i/ery  good  agreement. 

The  writei  is  deeply  indel  ted  to  Dr.  F.  L.  Martin 
(Professor  of  Meteorology)  for  his  suggestions  and  continue: 
help  throughout  the  L  ■••  ej  tigatipns  and  I:  ring  preparation 
o!  this  paper.   Special  credit  is  d  e  Professor  Martin  for 
his  large  share  in  developing  the  derivations  Ln  this  stud] . 
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LIST  01      ■  S    USED 
Symbol  ie  fin  it  ion 

u  .)i:.  ■    !  peed 

u  Friction  Velocity 

k  Von  Karmann  Constant 

Height  above  Lhe  Grouau 
Deacon  profile  number 


P 


Richar  so  i  Number 
Acceleration  of  Gravity 


© 


L 


? 


Potential  Temperature 


c£-  i]  i-i      /  Constant 


Mon  Li  -0    Likho  .    Scale 
I  h 


c.  Specific    Heat    oj    Air 


Density   of  Air 

H  .  Li  a  I  Flux  ol"  Heat  uy 

Turbulent  Diffusion 

y   Diffusivity  for  Heat 
Conduction 

Eddy  Diffusivity  for  Heat 
Momentum 

S  Normalized  Logarithmic 

Wina  Shear 


Rf  FIu;^  Richardson  Num 

V 


Constant  Edoy  Stress  of 
t  K  Surface  Layer 


Ratio  of  Convective  to 
V"  Mechanical  Energy  Source 

of  Turbulent  Energy 


Symi  -)I  .  ■  t  Lon 

2,   o  trie  Mean  of 

z 

R*  nsional   Heat 

£  (KH/K;)y 

RMS  Root  Mea  i  Squ 


v 


1.   Introduction. 

The  Deacon  profile  number  may  bt    define* 


^   „  *V_    2*^ 


"5  2     k 


(i) 


where 


u.,.  -  the  friction  velocity 

k   =   the  Voi:  Karmann  constant  and  is  approximately 
equal  to  O.-s- 

z        -     the  height  above  the  ground 

j3   =  the  Deacon  profile  number 

An  alternative  definition  of  p  which  Ls  equivalent  to 

that  in  equation  (1),  provided  p  and  u...  are  constant 

with  height  in  any  part  of  the  surface  layer,  is 


6«v 


P 


-2 


^2' 


dik 


(2) 


With  the  use  oi:  finite  differences  in  the  layer  1,  2,  i , 
the  last  definition  of  w5  implies  that  an  integral  mean 


value  is  , ;  .  en 


P 


(3) 


In  the  derivation  of  equation  (3),  the  successive  levels 

Zi,  *:v  ani  Zq  r,        to  e  '  ou   ed"  levels  so  tnat 
23/20   =  '<■    i-        -  2.   Iquation  (.>)  aiforos  a  means  of 
verjn.yin-3  Lhe  ■, ai.  .-.:     3  cor.hnit.cu  ,.y  oti.er  methods  to 
be  discussed  later, 


The  Rich.       i  .       y 


where 


g   -  the  acceleration 

(j  -  the  potential  temperat  ire 
layer 

u  =  the  wine  sp<  ei 
The  Richardson  auml  er   pi  ;ents  th<  rat  3    '.    :ui 
energy  produced  j  atmospheric         To  that  . 
by  mechanical  frictioi  , 

Mi  iin  and  Oi  ukho\  1  -,  54  I   ha     it  oduc 
profile  which  has  the  fo: 

*a        He  I  i 

in  place  of    the   equation      04JL/o«£  11  ^/\ k 2 

in  strictly  neutral  conditions. 

In  equation  (5)  0^-  is  a  constant  esti      by  Mouii 

Obukhov  to  i  e  0.6,  and  Lj  is  the  Mcnin-i 

length  defined  according  to  the  pair  of  equations 

ko,  H 
o 


7 


H  -  ~  k  oc  -^- 


& 


and  P  Lc  heat  and  density 

,u.  ,  respectively,  a     .  i  edd      usivit)     heat 

conduction,   [n  the  s  rface  boundary  layer,  u_,  ana  H  arc 

.  erall)  re  ;ar<  ed  as     tant  ..ic'u  height,  although  for 
H   this  condition  i   more  appropriate  above  i  meter 

!   riesi    ,    39  |  .   It  Ls  .veil  known  chat  cc,  at  ion  ( 
applies  ■).. ,   i  i  near  neutral  conditions.   However,  Ellison 

J  357  ]  has  recently       (te<  a  more  general  relation- 
ship I      q  ation  (9")   elow  \    ,   which  covers 


■  r 


ai  •.  v    surface  -layer  it  Li  ies  , 

En  de:      the  di  theoretical  relationship,  the 

"normalize!      tr  thmic  Ln   hear  S  has  been  used.   S  is 
defined 

5  -  k£  kit 

•W^s  (7) 

Usinj   !,  ]  and  Ri,  as  equal  Lons  (7),  (  )  and 

(4)  respectively  ,  one  can  derive  th<  equation 


zr   — 


K  o\ 


Mt)5  =  -(Rf)S 


(8) 


,   ..    is    the    ratio   ox    edd>         .  cus   vities    for   .cat 
a;.  rr», 

Rf    is    the   so-called    flux   UiCiiarv.soi.   .ui:....cr . 

In   d<     i        ig   equation    (8)     me   makes   use   of    the  lown 

def initio  i  c  ErictL  - ocity 

„*  I 

*  p 

J         3 


rty      ^  OIL- 
-rr'^     {-'       -  ls   the  constant  eouy  st] 

* '  -  d£- 

t  .  surface  Layex*,  and  D  i.:;  the  density  of  air. 
On  the  other  hand,  Ellison  11957  '.as  designed  an 
interpolation  ormula  for  the  .find  profile  which  fits 
observed  data  under  certain  Limiting  conditions  of  sta- 
bility.  His  suij^estcu  formula,  after-  some  transforma- 
tions, i;as  the  (form 


5*-  BJS   r  L 


where  V'  is  the  ratio  of  convective  to  mechanical 

energ>  sources  o^  turbulent  energy.  By  utilizing  equa- 
tions (7)  and  (J),   one  obtains 

-V  -! 


s*  p-xw=l}-x 


(10) 


For  small  values  of  z/L  or  oi  Ri,  one  can  easily  verify 
that  equation  (i.0)  gives  t  te  same  form  upon   inomial 

r  i 

expansion  as  that  oi   on  ]  equation  (j) 

L  J 

com  tests  of  numerous  wind  profiles  at  various 
micrometeoro logical  sites,  Panofsky,  Blackadar  and 
.  fehil   i  L960  I  conclucie<  Erom  equation  (10)  that 

the  best  fit.  The  signifi- 


if"--  i^v  ;i-v's =  I8*° gii 


cance   of   this   value    /ill    be   seen   in  Section  4. 


Lchai 
r<  nee    techniqu  taining   Ri  at   a 

Level   Zi    ■,   ..c-a^   the  midpoint   of    Layer   Zi    to   :.>     *.vas 

;gestec    b>     Lettau     I  see   pp,  328-329,    Lettau  aiui 
Davidson,    L957  I     ,   and   has    the    fo 


(11) 


where  z  am  .:...  are  "successive  dou  Led"  Levels,  and 
'■'--.,      i-  the  geometric  mean  Level  o.-:  z,  and  ;...   defined 
by  ?   -)"  V«;r  2-  s.V£  £     In  equation  (11),  t/.  and  u 
apply  to  the  top  oi  the  .aye;:,  and  C^  and  u.   apply  to 
fhe  lower  boundary,  while  Ri  is  considered  applicable  at 

the  geometric  me;         z        .   fote  Li;ac  Ri  is  negative 

i  ..  2 

in  unstable  conditions. 

3.   The  treatment  of  the  data. 

The  data  employ  -  ,  tre  selected  from  Table  8.1  and 
8,2  of  the  record  oi  Project     Lrie  Grass  Barad  et  al, 
1958  \    and  were  restricted  to  cases  of  Eree  convection. 
Priestley',;  criterion  for  Eree  convection,  »»  Ha1  :~  \  C<0'/:., 
has  been  used,  wher<  t  e  su   cript  1.5  refers  io  a 
measurement  of  Ri  at  1.5  meters.   From  the  wind  eata  of 
Project  Prairie  Grass,  44  case:  \i    free  convection  were 
selected,  and  from  these  a  composite  wim  ,     Le  ..as 
compute. .   The  mean  wind  speeds  at  each  doubled  Level  are 
listed  Ln   Table  1.   For  tliese  same  Eree-convective  cases, 
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and  isi  cases  of  free  convection,  H  has  he  average  value 
0  .f     wever,  M,   may  be  expi      in  the  Eorm 


H^Skaf^ 


K 


M 


t\ 


( 1 2 ) 


i./>  apply Lng  equations  (4)  awe.  (~0)  with  the  first  form  of 
H.  .  Thus  by  equation  (12)  with  '!,..  •-  0   ,  the  (RijKjj/Kw) 
relationship  for  cases  j;  Eree  convection  ... 


;-'H  v-  -JLi   - 


c  /  o  ^ — 


using  the  .due  k  =  0,     for  the  von  .Car.oa;:  constant. 

r 

In  using  equation  (1.3?,  V  =  ('h,.Xj  Pf  =  18.0  was 
aysu.net.  to  hold  exactly  at  the  level  j..:i  meters,  that  is, 
ere  ..i  =  -(  L'J  1 .   This  gives  the  fo  Lowing  results  at 


i.  =  1.41  meters : 


^h  _ 


Hence  i'or  tht 

other  layer  centers  ol  'laoie  1  ant.  the  (.li,    K,,,-  K.  ) 

relationship  cecomes 


( 


^H 


'i::e  ratio  K../K,  Ls  thus  obtained  as  the  solution  of  the 
cul  ic  v-q.ati.tui  (14)  at  al3  ol  ..  gi  o  tetric  mean  Levels. 
Table  2   shows  the  values  of  Ku/K^  computed  bj  equatxoa  (14) 


at  the  ubl<     rels   usii      coi  pu  t  i 

Ricl  ;i4  L  , 

Since  PanoJ    et  a   j  L960  »  obta  i<     o<  a, ;ree- 
ment,  on  the  average,  using  v  =  (Ku/Km)  A"*-  18.0,  this 
value  oi  \*     ..as  assumed  to  hold  exactly  in  this  stui 
at  the  leve]   .  •  meters.   The  iiain  justification  of  this 
assumption  i       this  I      is  clos   -   tl      letric 
mean  of  the  height  range  o      th<  wini  li  als  i         sti' 
jation.   It  turns  out  that  the  choice  V'  =   I,  .0  at  L.41 
meters,  results  In  values  of  V"  slightlj  jreater  than 
18.0  above  L.    leters,  but  somewhat  smaller  below  this 
Level.  Moreover,  when  these  V"  values  have  been  plotted 
versus  log  z      Eor  25  cm  <1  2  Cx  1600  cm,  an  average 
value  of  \.^    very  close  to  L8.C  results. 
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ogarithmic  erent  La  respect    to     z      leads    to 


1  6b    _    i 


_&«    . 

-    'O 

/  0  2 

'  -  -        - 

Fror.i  equation  (2)  the  sec  m        Lthin  the  u   ithes es 
i   recognized  to  be  —  K  ,   Also  upon  applying  logarithmic 


clif ferent .:.a<. ion  to  the  normalizec  lo  arithmic  wind  shear 

/ 

S   =   !  i  •"  ,v  •  .       one  obta  l  is 

i   *>5  _     L  ;  ■'•■  i: 

■4 

Elimination  of    -  ~   \      Erom  equations  (lii)  auu 
(16)  leads  to 

1       x   ■       K    v 


3 


assuming  u   t(     constant  Ln  the  layer. 

Hence  integration  of  equation  (17)  Erom  :. ,  Lo  ;•..,  yields 

the  mean  value  £j  for  the  layer  given  uy 
\ 

ft-  .  .{l  A'  --J.    / 


(17) 


-  ) 


where  K.j/lC,   is  also  a  function  of  Ri,   This  equation  is 

then  the  theoretical  p  —  ki  relationship.   The  value 

—     if — 

>'2l       is  to  e  associated  with  the  geometric  mean  height 

r 

ii 


z 


_   =  i  t    .  |  I  i!  .._;   Application 


of  eq  ati     v.;  ayers       n  successive 

doul  -       -  oi  la  >le  2  at  which  a  value  of  i',,,  •  r, 
is  !:r:own. 

Since  V"  =  L6.49  Ls  treated  as  a  constant  within 
the  surface  layer,  equation  (18)  becomes 

p     =    rt  u ;        __J_L_ _...>^  }■     (19) 

The  results  obtained  using  t::is  equation,  together  with 

the  t.ata  of  Table  i  and  2,  are  listei  in  Table  3.  In 

this  table,  the  foli   Lng  three  different  types  of  fo 

values  are  tabulate*  against  elevation: 

j3  i     ,  the  value  resulting  from  equation  (19) 
with  Ku/K^j  values  taken  from  Tc/ij   2. 

t  •  <-  ,    the  value  resulting  :rom  equation  (18) 

with(Ku/Kj,j)  q  =  18,0  at  all  levels,  so  that 
satisfies 

6-i  -^O-'8  *     '  'l1'"      ■■     (20) 

- 

,  the  value  resulting  Erom  equation  (3). 

pss  ~ 

It  should  be  noted  that  G   an<   pT   are  in  reia- 
tively  close  agreement,  but  some  unusual  discrepancies 
occL.r  in   /  x    .   This  difficulty  is  dealt  with  at  more 


Length  in  Section  6. 
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6.   Correct:     c  /3    values. 

Some  inconsistency  occurs  in  computing  fa  using 
direct  wind  <a...a.   At  the  200  cm  level,  the  ft    value 
was  very  low  in  comparison  with  the  surrounding  |3  's ■ 
This  error  occurred  because  the  vertical  .;ind  spee 
increment  between  the  200  em  and  400  cm  level  was  not 
consistent  with  those  in  adjacent  layers.  This  suggests 
that  there  may  have  been  an  instrumental  error  at  one 
or  more  ol"  these  levels.   Hence  a  control  '.a  ta  tai-  le  was 
made  up  by  selecting  :;ata  occurring  at  the  time  oi  a 
neutral  wind  profile,  that  is,  when  the  potential  tempera- 
ture is  isothermal,  with  height.  This  usually  occurs  near 
sunrise  and  sunset.  Fourteen  (14)  cases  have  been  selected 
that  meet  those  requirements  and  from  these  a  near -neutral 
wind  profile  was  computed,  as  shown  in  Table  4. 

A  perfect  neutral  wind  profile  would  have  the  char- 
acteristic of  constant  (u2  -  ui)  increments  between 
successive  dou  led  levels.   However,  lI;c  values  of 
(u2  -  ui)  in  Table  k    indicates  that  this  was  not  the 
case.  The  lowest  three  layers  in  Table  4  have  (uo  -  u^) 
increments  which  are  quite  close  to  their  overall  mean 
of  49.6  cm  sec   »   This  last  value  indicates  the  -lope 
of  the  logarithmic  profile  which  exist:.;  under  neutral 
conditions.   The  remaining  three  wind  increment:;  of 
Table  4  then  indicate  percentagewise  how  much  the  next 
three  layers  deviated  from  consistency  with  the  lower 
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,  ,  a      ct ion  factor 
has  ,    each  Layer.   For  example,  th  cor- 

rection Eactoi  0.  i2i  Ls  bas    m  the  fact  that  tht  nean 
wind  increment  .  ico  ciei  Ls  60.3  cm  see  "  rather  than 
49.6  cm  sec   .   In  arriving  at  the  corrected  wind-speed 
differences  of  Table  y,    on<  must  enter  with  the  data 
Tables  L  an;      Vhir.  cor  example,  the  mean  ./ind      rence 
in  the  layer  200  to  400  cm  Ls  i   idicated  ay  48.1  cm  sec 
in  the  non-ne  tral  case  of  Table  1.   However,  the  data  of 
Table  4  Indicates  that  this  reading  is  overestimated  in 
the  neutral  case  _   he  factor  1/0.822.   Hence,  the  cor- 
rected wind  speed  increment  applicable  to  the  same  layer 

is 

-i 
0.02Z   .:   48.1    cm   sec        =  39.5/   cm  sec 

Similar   corrections   ma>    then     e  applie<     to    the   other 

layers . 

With  the  corrected  wind  differences  o. taiaeu  bj  this 

method  and   Lspla>      Tal  Le  5,  correctei   n_  val 

are  then  obtain*        .,  at  ion  (.3;.  and  are  Listed  in 

the  last  ro  "  of  the  ta  >le. 
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Finally  as  a  c(      on,  c  imparison 
/a lues  obtained  t 

directlj  from  the  correcte<  .vind  data  Indicates  tl 
tin.1  ....    tical  ": ..   'ell  subst an t Latec . 

The  average  ciiff  srence    tween   |3T  and 
is  eq;;ai  to   -0.031  with  a^  ,..  S  erroi  of  0.094. 
Similarly,  the  average   LJ  ere  ce  i  i  tv/een  I  -5~--  and 
H_  is    equal  1:j   -0.  >46  with  an  RMS  error  of 

Of  t  e  two  theoretical  formulas,  there  seems  to 
be  a  slight  advantage  .u  U.e  w...  of  equal  Lon  (20)  over 
equation  (19),  except  possible  eiou  25  cm,  where 
t  i  Lack  of  data,  comparison  is    not  possible. 
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